A Simple Cellular Automaton Model for Influenza A Viral Infections 



Catherine Beauchemin 
Department of Physics, University of Alberta, Edmonton, AB, T6G 2J1, Canada* 

John Samuel 
Faculty of Pharmacy and Pharmaceutical Sciences, 
University of Alberta, Edmonton, AB, T6G 2N8, Canada 

Jack Tuszynski 

. . . Department of Physics, University of Alberta, Edmonton, AB, T6G 2.J1, Canada 

'nJ- I (Dated: February 6, 2008) 

, Viral kinetics have been extensively studied in the past through the use of spatially homogeneous 

^SJ . ordinary differential equations describing the time evolution of the diseased state. However, spatial 

characteristics such as localized populations of dead cells might adversely affect the spread of infec- 
tion, similar to the manner in which a counter-fire can stop a forest fire from spreading. In order to 
investigate the infiuence of spatial heterogeneities on viral spread, a simple 2-D cellular automaton 
(CA) model of a viral infection has been developed. In this initial phase of the investigation, the 
' CA model is validated against clinical immunological data for uncomplicated infiuenza A infections. 

, Our results will be shown and discussed. 

; I. INTRODUCTION 

d : 

• ^ In the past two decades, many approaches have been chosen to model different aspects of the immune system. 

iP ■ Differential equation models are perhaps the most common and are typically used to simulate the immunological and 
O^l epidemiological dynamics of particular diseases to try and identify the critical parameters involved 0, . Genetic 

■ algorithms have been applied to the modelling of the evolution of diversity and pattern recognition capability in the 
I \ immune system @,^3E3- On the theoretical front, some ingenious work j^ [T5li20,. .22] has been done to understand 

^ . aspects of the immune system from the perspective of optimization problems. Efforts to build an immune system 
fS| ' tailored for computer networks l2ll| have raised our understanding of the crucial nature of certain immune system 
T-H , mechanisms. Cellular automaton (CA) models have been developed to simulate the immune dynamics of particular 

■ diseases jl^l, to model shape space interactions |2^ based on the network theory introduced by Jerne and have 
] been chosen as the tool of choice to implement highly complex general immune system simulators 0, l4l Il4l ITgj that 

are to be used to run immune system experiments in machina. More recently, a promising stage-structured modelling 
approach was used to study cytotoxic T lymphocyte response to antigen 

Of all of the above, ordinary differential equations are the most commonly used method to model disease dynamics. 
, y By using ordinary differential equations — which describe a system assuming it is homogeneous in space — to study 
viral kinetic, one chooses not to take into consideration spatial characteristics that could potentially play a nontrivial 
' I role in the development and outcome of a viral infection. Here, the authors introduce a cellular automaton model 
. . that will be used in later work to study the influence of spatial heterogeneities on the dynamical evolution of a viral 
^ infection. In this first step of the project, the CA model is used to simulate an uncomplicated infection with the 
k>( influenza A virus in order to verify that it is not only capable of replicating the typical shape of an immune response 
j_] to an uncomplicated viral infection, but that it can also give quantitatively reasonable results when parameterized 
' for a particular viral infection. In Section m we give a brief description of the nature of influenza. In Section IlIII we 
describe the structure and evolution rules of the CA model. In Sect ion HVl we present results from the CA simulations 
and compare them against data from the literature. Finally, in Section^ we expand on the biological meaning of the 
model's parameters, how we have arrived at the values we chose for the model, and examine the model's sensitivity 
to their values. 



o 
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II. ABOUT INFLUENZA A 



Influenza, in humans, is caused by a virus that attacks mainly the upper respiratory tract, the nose, throat and 
bronchi and rarely also the lungs. According to the World Health Organization (WHO), the annual influenza epidemics 
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affect from 5% to 15% of the population and are thought to result in between three and five million cases of severe 
illness and between 250 000 and 500 000 deaths every year around the world. Most deaths currently associated with 
influenza in industrialized countries occur among the elderly over 65 years of age |25| . 

Influenza viruses are divided into three groups: A, B, and C. Additionally, influenza viruses are defined by two 
different protein components, known as antigens, on the surface of the virus. They are spike-like features called 
haemagglutinin (H) and neuraminidase (N) components. Influenza A has two subtypes which are important for 
humans: A(H3N2) and A(HINI), of which the former is currently associated with most deaths. Antibodies to H 
are strain-specific and neutralize the infectivity of the influenza A virus, while antibodies to N have a less protective 
effect. Thus, it is suggested in 3J that the the antigenic properties of influenza A virus be associated with the H 
determinant. 

It is the high mutation capability of the influenza virus that makes it a great public health concern. The genetic 
makeup of influenza viruses allows for frequent minor genetic mutations to take place. This makes it necessary to 
constantly monitor the global influenza situation in order to adjust the influenza vaccines' virus composition annually 
to include the most recent circulating influenza A(H3N2), A(HINI) and influenza B viruses. Additionally, influenza 
A viruses, including subtypes from different species, can swap, reassort, and merge genetic material resulting in novel 
subtypes. Three times in the last century, influenza A viruses have undergone major genetic changes mainly in their 
H-component, resulting in global pandemics and large tolls in terms of both disease and deaths |24|. 



III. THE CELLULAR AUTOMATON MODEL 



Our CA model considers two species: epithelial cells, which are the target of the viral infection, and immune cells 
which fight the infection. The virus particles themselves are not explicitly considered, rather the infection is modelled 
as spreading directly from one epithelial cell to another. The CA is run on a two-dimensional square lattice where 
each site represents one epithelial cell. Immune cells are mobile, moving from one lattice site to another, and their 
population size is not constant. The CA lattice is therefore like a tissue of immobile cells which is patrolled by the 
mobile immune cells. The CA is updated synchronously (i.e. all cells are time-stepped at once rather than one cell at 
a time). The boundary conditions for both the epithelial and immune cells are toroidal, i.e. an immune cell moving off 
one edge of the grid are reintroduced at the opposite edge and an infectious epithelial cell at one edge of the grid can 
infect healthy cells located at the opposite edge. Finally, the neighbourhood of a lattice site is defined as consisting of 
the 8 closest sites. Infected epithelial cells can only infect their 8 surrounding neighbours and immune cells can move 
at a speed of 1 site per time step to any of the 8 neighbour sites or remain in place. Below are the detailed evolution 
rules for each of the cell species. 

Epithelial Cells An epithelial cell can be in any of five states: healthy, infected, expressing, infectious, or dead. 
Transitions between epithelial cell states occur as follows: 

• Epithelial cells of all states become dead when they are older than CELL_LIFESPAN. 

• A healthy epithelial cell becomes infected with probability INFECT_RATE/(8 nearest neighbours) for each 
infectious nearest neighbour. 

• An infected cell becomes expressing, i.e. begins expressing the viral peptide, after having been infected for 
EXPRESS .DELAY. 

• An expressing cell becomes infectious after having been infected for INFECTJ)ELAY > EXPRESSJDELAY. 

• Infected, expressing, and infectious cells become dead after having been infected for INFECTXIFESPAN. 

• Expressing and infectious cells become dead when "recognized" by an immune cell. See below for the 
meaning of "recognition." 

• A dead cell is revived at a rate DIVISION_TIME~^ x # healthy/^ dead. When revived, a dead cell becomes 
a healthy cell or an infected cell with probability INFECTJlATE/(8 nearest neighbours) for each infectious 
nearest neighbour. 

A simulation is initialized with each epithelial cell being assigned a random age between and CELLXIFESPAN 
inclusively. All epithelial cells start in the healthy state with the exception of a fraction INFECT_INIT of the 
total number of epithelial cells which, chosen at random, arc set to the infected state. 

Immune Cells An immune cell can be in any of two states: virgin or mature. A virgin cell is an immune cell that 
has no specificity. A mature cell is an immune cell that has either already encountered an infected cell or has 
been recruited by another mature immune cell. 
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• Immune cells move randomly on the CA lattice at a speed of one lattice site per time step. 

• An immune cell is removed if it is older than IMM_LIFESPAN. 

• An encounter between an immune cell and an expressing or infectious epithelial cell requires the immune 
cell to be in the same site as the expressing or infectious cell. 

• A virgin immune cell becomes a mature immune cell if the lattice site it is occupying is in the expressing 
or infectious states. 

• A mature immune cell occupying an expressing or infectious lattice site "recognizes" the epithelial cell and 
causes it to become dead. 

• Each recognition event causes RECRUITMENT mature immune cells to be added at random sites on the CA 
lattice after a delay of RECRUIT_DELAY. 

• Virgin immune cells are added at random lattice sites as needed to maintain a minimum density of 
BASE_IMM_CELL virgin immune cells. 

A simulation is initialized with a density of BASE_IMM_CELL virgin immune cells at random locations on the CA 
lattice, each with a random age. 

Implicit in this model are the following physiological assumptions: 

• Only healthy epithelial cells are able to divide. 

• Immune cells cannot be infected and so only ever die of old age. 

• No memory immune cells are considered. Memory could be added by setting an extended lifespan to a given 
proportion of immune cells created during a viral invasion. 

• We consider a single viral strain with no mutation. This implies the existence of a single, unique, epitope 
identifying the infection. 



IV. RESULTS 



A great deal is known about the influenza virus including viral structure and composition, the replication process, 
and even some dynamical data regarding the viral and antibody titer over the course of the infection P, |^, [l3l |. 
However, key dynamical information such as the flow rate of immune cells within an infected tertiary lymphoid organ, 
the clearance rate of viral particles, the lifespan of an infected epithelial cell, is either uncorroborated, unknown, or 
known with poor precision. In our model, all but two parameters have been taken directly or adapted from Q. They 
are presented in Tabled along with their value and description. 

Due to the limited dynamical information available for influenza A, there are only a few quantitative characteristics 
against which we can compare our model. They are: 

1. The infection should peak on day 2 (48 h) 

2. Over the course of the infection, the fraction of epithelial cells that are dead should be as follows 

(a) 10% on day 1 (24 h); 

(b) 40% on day 2 (48 h); 

(c) 10% on day 5 (120 h). 

3. From ^3], virus concentration should decline to inoculation level on day 6 (144 h). From experimental data 
recovered from 8 volunteers indicated that virus shedding persisted for 5 ± 2 d (72-168 h). 

4. The number of immune cells should peak anywhere between day 2 (macrophages' peak) and day 7 (cytotoxic T 
cells' and B cells' peak) (48 h-168 h) ^. 

5. At their peak, the number of B cells, helper T cells, and cytotoxic T cells should be 100-fold greater than their 
normal concentration, while that of plasma cells should be lO^'-fold greater Q. 
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Parameter 


Value 


Description 


gr i d_w i d"t h. 


440 


Width of the grid (ceUs) 


gir i d._li6 i ght 


280 


Height of the grid (cehs) 


FLOW_RATE 


6 ts/h 


Speed of immune cells (time steps/hour) 


TMM T TFFSPAN 


168 h 


T,ifpt;nflTi nf an iTTiTTiiiTiP ppII 


CELL_LIFESPAN 


380 h 


Lifespan of a healthy epithelial cell 


INFECT_LIFESPAN 


24 h 


Lifespan of an infected epithelial cell 


TNFFrT TMTT 


0.01 


n V a pi" 1 riTi Or inii"ifill'\7' iTiTpr*i"pri ppIIq 

± Idl^LjiWll KJl. llllLifXliy llllCVyLiCU- VjCliO 


INFECT_RATE 


2h-i 


Rate of infection of neighbours 


INFECT_DELAY 


6h 


Delay from infected to infectious 


EXPRESSJ3ELAY 


4h 


Delay from infected to viral expression 


DIVISION.TIME 


12 h 


Duration of an epithelial cell division 



BASE_IMM_CELL 1.5 x 10~'' Minimum density of immune cells per epithelial cell 
RECRUIT_DELAY 7 h Delay between the recruitment call and the addition of immune cells 
RECRUITMENT 0.25 Number of immune cells recruited when one recognizes the virus 



TABLE I: The model's default parameters extracted from Q. 
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FIG. 1: Simulation results averaged over 10 simulation runs using the parameter set presented in Table Q] The paired lines 
mark one standard deviation and represent the proportion of epithelial cells that are healthy (dashes) , containing the influenza 
A virus particles (dots), dead (solid), as well as the proportion of immune cells per epithelial cell (dots and dashes). 

The simulation results for the default parameters presented in Table are shown in Figure ^ The infection peaks 
on day 2, and although there are fewer than 10% of the cells dead on days 1 and 5, there is 40% of the cells dead 
on day 2. Recovery from the infection is faster than experimental data suggests, but this is probably mainly due to 
the fact that immune cells in our model are generic. During a true influenza A infection, the first part of the immune 
response is dominated by macrophages and T cells while the later part is dominated by the action of antibodies, a 
subtelty that is missed by our bulk model. The number of immune cell does peak between day 2 and 7, with an 
approximate 1000-fold increase from the normal concentration set by BASE_IMM_CELL. 

V. DISCUSSION OF THE PARAMETERS OF THE MODEL 

Some of the model's parameters have quite a large uncertainty. In some cases this is a reflection of the difficulty of 
measuring the parameter while in other cases it is a reflection of the variation between individuals. The model also 
has two free parameters, RECRUITMENT and INFECT_RATE, whose values are not known from physiological data. Here, 
we will discuss each parameter in detail and investigate the sensitivity of the model their variation. 
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FIG. 2: Effects of the lattice size on the dynamics of the infection. The greyed areas represent one standard deviation after 
10 runs while increasing darkness represents increasing grid sizes of 7,700 cells, 30,800 cells, 123,200 cells, 1,971,200 cells, and 
9,432,500 cells. For clarity, the inset shows a detail of the healthy cell curves alone. Increasing the grid size is seen to not affect 
the average behaviour but decreases the standard deviation. 

A. grid_width and grid_height 

Since each lattice site represents an epithelial cell and the width and height of the grid are given in numbers 
of epithelial cells, it is possible to calculate the physical size of the simulation. The model's default lattice has 
grid_width x gridJieight — 123, 200 cell sites. In uninflamed tissue, there are about 2,200 ciliated cells per 
square millimetre of epithelial area Q. Taking ciliated cells — which account for 60% ~ 80% of internal bronchial 
surfaces — to correspond to each lattice site of the CA model, the real-life size of the default lattice is about 
123200 cells/2200 cells • mm~^ = 56 mm^. When inflamed, there are about 700 cells per square millimetre of epithe- 
lial tissue area, giving the simulation a real-life area of 176 mm^. 

To test the sensitivity of the model's algorithm to changes in lattice size, i.e. to determine if the finite nature of 
the simulation influences the results, simulations have been run with lattices ranging from 7,700 cells to 9,432,500 
cells (the most that could be simulated with the 512 MiB of RAM available in the simulation machine). The results 
obtained for several lattice sizes in this range are presented in Figure [3 What one sees from this plot is that the 
average behaviour of the model is unchanged by modifications to the size of the simulation area, but that the standard 
deviation of the behaviour about the mean decreases as the simulation area is increased. This is consistent with the 
interpretation of the results from a single simulation run on a large grid as being equivalent to averaged results from 
several simulations run on smaller grids. This would be the case — that a large simulation is equivalent to an average 
over several smaller simulations — anytime the simulation size itself does not play a role in the dynamics, and so 
observing this behaviour strongly suggests that the default lattice size is sufficiently large for its finite nature to be 
unimportant. 

B. FLOW_RATE 

Due to the discrete nature of a cellular automaton, immune cells move by increments of one lattice site per time 
step. Because each lattice site represents one epithelial cell, each time step has to represent the true time required 
for an immune cell to go from one epithelial cell to the next. Very little is known about the speed at which immune 
cells patrol tertiary lymphoid organs or by how much the flow will vary depending on whether the area is healthy or 
infected. This has made the model's time scale difficult to establish. 

All time dependent parameters in the model are scaled by a parameter named FLOW_RATE which is defined as the 
number of CA time steps corresponding to 1 h in real time. The duration of the CTL lethal hit process, i.e. the time 
it takes for a cytotoxic T lymphocyte to kill an infected cell, is about ^attack = 20 min ~ 30 min Q . One can suppose 
that an immune cell that is not involved with, or in the process of destroying, an infected cell will move from one 
epithelial cell to the next in less time, say ifrcc- Since on average about 15% of cells are infected over the course of 
the infection for the default set of parameters chosen for the model (see Figure 0, one can compute the average time 
it takes for an immune cell to go from one epithelial cell to the next as tavg = 0.15 x iattack + 0.85 x thcc- Assuming 
that tfrco can be anywhere between negligible, « min, to at most 30 min, then iavg can range from 3 min to 30 min, 
which corresponds to a range of 2 to 20 time steps per hour for FLOW_RATE. We have chosen to set FLOW_RATE to 
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FIG. 3: The effect of varying FLOW_RATE on the viral infection's dynamics. From left to right, the graphs show the behaviour 
obtained using FLOW_RATE values of 2, 6, and 20 time steps/h respectively (the physiologically plausible range). The central 
graph is the same as Figure Q These show that the simulation is sensitive to the parameter FLOW_RATE. 




192 216 240 264 



FIG. 4: The effect of varying IMM_LIFESPAN on the viral infection's dynamics. From left to right, the graphs show the behaviour 
obtained using IMM_LIFESPAN values of 48 h, 168 h, and 480 h respectively. The central graph is the same as Figure These 
show that the simulation is only sensitive to the parameter IMM_LIFESPAN when its value is smaller than the time required for 
the immune cells to clear the infection. 



6 time steps/h, which means that immune cells move at a speed of one epithelial cell per 10 min. Note that this 
corresponds to a tf^cc of about 6 min ^ 8 min, but we have found no experimental value in the literature to compare 
this to. 

Figure 131 shows the effect of varying FLOW_RATE throughout this range. The faster the immune cells move (the larger 
FLOW_RATE is), the bigger the advantage they will have over the growth of the infection and thus the less pronounced 
the infection. This is clearly illustrated in Figure O where the smallest value of FLOW_RATE results in a lethal infection 
while the largest value results in a minor infection. The CA model is extremely sensitive to this particular parameter 
and therefore it would be important for future work to determine this parameter with greater precision. 



C. IMM_LIFESPAN 



The problem in defining this particular parameter with much precision comes from the fact that the immune cells 
in our CA model are generic, i.e. they do not represent a particular species of immune cell such as cytotoxic T cells 
or plasma B cells, rather they represent the combined action of all involved species of immune cells. In 01, the rate 
constants of natural death for immune cells range from 0.05 for B cells to 0.5 d^^ for cytotoxic T cells and plasma 
B cells. This corresponds to an immune lifespan range of 48 h ~ 480 h. For our model, we chose an immune cell 
lifespan of 168 h as the default value. In Figure 01 the dynamics resulting from the use of a 48 h and 480 h immune 
cell lifespan is shown. In broadest terms there are two outcomes and which one is observed depends on whether the 
immune cell lifespan is greater than or less than the time required for the immune cells to completely eliminate the 
infection. When the immune cell lifespan is less than the time required to clear the infection, all cell populations 
oscillate, much in the manner of a predator-prey system — a situation which may be referred to as a chronic infection. 
When the immune cell lifespan is greater than the time required to clear the infection, it does just that. 

From Figure 0] it can be seen that only at the very bottom end of the physiologically plausible range for the 
parameter IMM_LIFESPAN does its value begin to influence the model's dynamics. Over most of this parameter's 
range, the infection is killed too early for the immune cell lifespan to be important. Since experimental data indicates 



FIG. 5: The effect of varying DIVISION_TIME on the viral infection's dynamics. From left to right, the graphs show the behaviour 
obtained using DIVISION_TIME values of 7 h, 12 h, and 24 h respectively. The central graph is the same as Figure These 
show that the simulation is sensitive to the parameter DIVISION_TIME. 



that the virion concentration should have declined to inoculation level no later than 144 h (6 d) into the infection, we 
chose a value for the immune lifespan that is slightly larger than this, namely 168 h (7 d). 



D. DIVISION.TIME 

The duration of a single division of an epithelial cell is about 0.3 d ~ 1 d Q- We chose a value of 12 h (0.5 d) for the 
model parameter DIVISIONJTIME. Assuming that only healthy cells can undergo successful division, the probability 
per unit time that any dead cell is revived is given by 

r,/ , , ,. X 1 # healthy cells 

-Pidead alive) = x — — - — — . 

^ ^ DIVISIDN.TIME # dead cells 

Figure El shows the effect of varying DIVISIDN_TIME over its physiologically allowed range of 7 h to 24 h. One 
sees that this parameter can have a significant effect on the dynamics and even on the outcome of the infection. 
DIVISION_TIME represents the capacity of the system to regenerate itself: the greater the regenerative power of the 
system — the smaller DIVISION_TIME — the better its chances of recuperating after a severe infection. In this model, 
at the high end of DIVISION_TIME's range there is a tendency for the influenza A infection to kill the host, while at 
the low end of its range the epithelial cell regeneration rate is too high to reproduce the level of damage observed 
experimentally. We chose a default value of 12 h for our model. With this value, 40% of the epithelial cells are dead 
at the peak of the infection which is in agreement with the experimental data. 



E. CELL_LIFESPAN 



An epithelial cell, in vitro, can go through 20 to 25 population doublings in its lifespan jl8| . Each population 
doubling lasts for DIVISIDN_TIME, so the hfespan of an cpithchal cell should be about 160 h ~ 600 h. We chose 
the mid-range value of 380 h as the default for our model. One impediment to finding more precise information 
about epithelial cell lifespans is that most in vitro research is done using cell lines that have been made immortal. 
Additionally, there is much variability in cell lifespan across cell types. In any case, even the lower bound of 160 h 
is longer than the duration of an uncomplicated influenza A infection and thus the value of CELL_LIFESPAN, within 
the physiologically plausible range, should have little effect on the results. Figure shows the dynamics of the viral 
infection for cell lifespans at the physiologically allowed extrema of 160 h and 600 h as well as at the default mid- 
range value of 380 h. In the model, one consequence of the en masse killing of epithelial cells by the viral infection 
is the synchronization of the ages of newly-grown epithelial cells. This is seen as a cell-death "echo" approximately 
CELL_LIFESPAN h after the peak of the infection. The only visible difference between simulations run with different 
values of CELL_LIFESPAN is a change in the time at which the cell-death echo occurs. 



F. EXPRESS_DELAY, INFECT_DELAY, and INFECT.LIFESPAN 



When an epithelial cell becomes infected by an influenza A virion, 4 h (EXPRESS_DELAY) elapse before it can start 
expressing viral antigen and thus be recognized by an immune cell . It will take the infected cell an additional 2 h 



8 




FIG. 6: The effect of varying CELL_LIFESPAN on the viral infection's dynamics. From left to right, the graphs show the behaviour 
obtained using CELL_LIFESPAN values of 160 h, 380 h, and 600 h respectively. The central graph is the same as Figure^ These 
show that the simulation is not sensitive to the parameter CELL_LIFESPAN. 
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FIG. 7; The effect of varying INFECT_INIT on the viral infection's dynamics. From left to right, the graphs show the behaviour 
obtained using INFECT_INIT values of 0.001, 0.01, and 0.1 respectively. The central graph is the same as Figure^ These show 
that the simulation is sensitive to the parameter INFECT_INIT. 



(INFECT_DELAY = 6 h) to begin releasing virus particles and thus become infectious Ultimately, influenza A virus 
being highly cytopathic, the Hfespan of an infected epithelial cell is shortened to at most 1 day (INFECT_LIFESPAN = 
24 h) 0. 



G. INFECT.INIT 



This parameter represents the magnitude of the initial dose of influenza A virus delivered to test volunteers in 
viral replication experiments. In particular, in our model this is the fraction of epithelial cells initially set in the 
infected state. This parameter's value is not bounded physiologically, rather it is an arbitrary initial condition. In 
studies of the influenza A virus reported in Q , the typical initial concentration of aerosol-delivered virions is about 
10^ virions/mL. Since the average concentration of epithelial cells in the area typically affected by influenza varies 
from 10^ cells/mL to 10^*^ cells/mL, and a single epithelial cell can absorb from 1 to 10 influenza virions, then the 
typical fraction of cells initially infected is computed to be 

10^ virions/mL 
INFECT_INIT ' 



(109 - IQio) cells/mL x (1 - 10) virions/cell 
0.001 - 0.1. 



We have chosen a default value of 0.01 — the logarithmic middle of the range. Figure presents the efl'ect on 
the model's dynamics of varying INFECT_INIT over this range. One sees that the evolution and the outcome of the 
infection is dependent upon the initial infective dose received. This reflects the competitive nature of the disease and 
immune system interaction: if one starts with too large a numerical advantage over the other, the latter is bound to 
loose. 
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FIG. 8: The effect of varying INFECT_RATE on the viral infection's dynamics. From left to right, the graphs show the behaviour 
obtained using INFECT_RATE values of 1, 2, and 4 infected neighbour(s) • h~^ respectively. The central graph is the same as 
Figure Q These show that the simulation is sensitive to the parameter INFECT_RATE. 

H. INFECT_RATE 

The growth rate of the infection in our model is determined by the parameter INFECT_RATE, which is defined as the 
number of its neighbours that an infectious cell will infect per hour. This parameter was not found in the literature as 
it does not correspond directly to the true biological process of infection. In fact, 6 h after a cell has been infected by 
the influenza A virus, it begins to release influenza A virions (virus particles) into its surrounding and it is those virions 
that move on to infect neighbouring cells. In our model, we have decided not to consider virions at all and instead have 
modelled the infection process by the direct infection of neighbouring cells by infectious cells. This is computationally 
simpler and more efficient. Throughout influenza A literature, however, the infection growth is described in terms 
of the viral titer over the course of the infection rather then as a growth rate. This makes it difficult to establish 
the parameter INFECT_RATE from experimental data. We have chosen to make it a free parameter and found that we 
obtained the best results by setting it to 2.0 h~^, which means that an infectious cell will, on average, infect 2 of its 8 
nearest neighbours per hour. One expects that a smaller value of this parameter will yield a less pronounced infection 
while larger values will increase the amplitude of the infection and potentially result in the death of the host. This 
is illustrated in Figure |S1 where we have explored the effect of using a value 1/2 x (1 infected neighbour • h^^) and a 
value 2x (4 infected neighbours • h^^) the default. 

I. BASE_IMM_CELL 

The viral infection model only concerns itself with immune cells that have a role to play in the viral infection, i.e. 
those immune cells with a receptor capable of recognizing an epitope of the viral strain present in the simulation. 
There must always remain a minimum density of cells capable of recognizing and responding to a new outbreak, 
and this density is set by the parameter BASE_IMM_CELL. The value for BASE_IMM_CELL was calculated based on the 
following information. It has been determined |2^J that intraepithelial lymphocytes have a frequency of about 15 per 
100 epithelial cells. Further, the probability that an immune cell chosen at random will recognize an epitope chosen 
at random is typically 10~^ [T^ . Finally, taking into account memory and cross-reactivity phenomena in the case of 
influenza A in adults, it is suggested in ^3] that the fraction of B or T cells capable of responding to a viral infection 
with influenza A be increased by 100-fold. Thus, BASE_IMM_CELL is calculated to be (15/100) x 10"^ x 100 = 1.5 x lO"'' 
virgin immune cells per epithelial cell. 

J. RECRUITMENT 

When a mature immune cell encounters an infected epithelial cell presenting a viral peptide, the mature immune 
cell will "recruit" a number of new mature immune cells. The number of newly created mature immune cells is set 
through the parameter RECRUITMENT. This parameter is one of the two free parameters of the model, and it has no 
direct biological basis because the immune cells in our model are generic. This parameter represents the efficiency 
with which the immune system responds to the viral infection. We explored a range of values for RECRUITMENT, and 
found that a value of 0.25, i.e. a 25% chance of recruiting a mature immune cell upon each successful recognition, 
gives satisfying results. 



10 




time (h) lime (hi lime (h) 



FIG. 9: The effect of varying RECRUITMENT on the viral infection's dynamics. From left to right, the graphs show the behaviour 
obtained using RECRUITMENT values of 0.05, 0.25, and 1.25 respectively. The central graph is the same as Figure^ These show 
that the simulation is sensitive to the parameter RECRUITMENT. 




FIG. 10: The effect of varying RECRUIT_DELAY on the viral infection's dynamics. From left to right, the graphs show the 
behaviour obtained using RECRUIT_DELAY values of 2 h, 7 h, and 12 h respectively. The central graph is the same as Figure 
These show that the simulation is not sensitive to the parameter RECRUIT_DELAY. 

The typical immune cell division time is about 12 h ^ 24 h ,3], which corresponds to a recruitment ranging from 0.08 
to 0.16 immune cells per hour. Considering that an immune cell encounters 6 epitheUal cells per hour (FLOW_RATE — 6 
time steps per hour) and that on average about 15% of those cells are infected, then our model's effective recruitment 
rate is about 0.25 x 6 x 0.15 = 0.23 immune cells per hour, which is higher than the biological range described above. 
However, that recruitment range is based soly on the cellular divisions of activated B and T cells and does not take 
into consideration the high production of antibodies by plasma B cells which contributes to reducing the infection, 
a factor that can account for the higher recruitment rate value required by our model. Simulation runs using a 
RECRUITMENT value of 0.05 and 1.25 are presented in FigurelHI One can see that the parameter RECRUITMENT has an 
important effect on the development and the outcome of the infection: too small a value makes the infection fatal, 
while larger values decrease the amplitude and duration of the infection. It is thus a key parameter and much effort 
should be invested in finding a way to draw a parallel between it and true biological parameters. 

K. RECRUIT_DELAY 

Once a "recruitment call" has been made, immune cells are not added immediately to the simulation lattice. The 
delay between the call and the actual addition of immune cells corresponds roughly to the typical transfer time of 
cytotoxic T lymphocytes from their main area of production — lung associated lymphoid tissue (LALT) compartment 
— to the affected area — the mucous compartment — which, in ^3], is said to take from 2 h to 12 h. We chose to use 
the mid-range value of 7 h RECRUITJDELAY's default value. Figure [TOl shows the effect of varying RECRUIT_DELAY over 
the biologically acceptable range. It is clear that changes to the value of this parameter have little influence on the 
dynamics of the model over the uncertainty range. 
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VI. CONCLUSION 

Here, we have introduced a CA model for an uncomplicated viral infection. We have shown that once parameterized 
for the particular case of influenza A, our CA model — which is described by 7 state variables and 12 parameters — 
is sophisticated enough to reproduce the basic dynamical features of the infection. One expects a 12 parameter model 
to be able to match 7 dynamical features, however all but 5 parameters of this model are sufficiently-well bound by 
physiological data that they cannot be used to tune its behaviour. Remarkably, our 7 variable and 12 parameter 
model's agreement with the experimental dynamical characteristics compares well to that of the 13 variable and 60 
parameter model presented in To the authors' knowledge, the CA model presented here along with the ODE 
model presented in ^ are the only two existing immunological models for influenza A viral infections. This is 
quite surprising considering its significance to human health worldwide, particularly its high potential for causing 
devastating epidemics j25j |. 

The results obtained with the model described herein have surpassed our expectations. We had hoped to develop 
a model that would be good enough to be used as a test bench to investigate various theoretical aspects of viral 
infections, and we have obtained an influenza A model that behaves impressively well when compared against available 
experimental data. We believe that by adding additional details to the model such as specific immune cell types, our 
CA model could be a very promising model of influenza A. 

A wide range of behaviours were observed as the model's parameters were varied over their biologically allowed 
range. It would be valuable to test whether these predicted behaviours can actually be observed in vivo, since if they 
are not seen then that would falsify the model. 

In future work, we will make use of the CA model introduced here to investigate various theoretical aspects of 
viral infections, for example the extent to which such things as spatial inhomogeneities and age classes can affect the 
evolution and outcome of a viral infection. It will be interesting to compare this CA model against the recent work in 
Pi which makes use of a stochastic stage-structured model, an alternative to both differential equations and cellular 
automata. 
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